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POST STALL DYNAMICS OF AIRCRAFT can produce 
uncontrollable motions which are often accom- 
panied by high yawing and rolling rates. If 
the longitudinal pitching moment is nose up, 
these initial motions are called incipient 
spins. If the pitching moment is nose down, 
they are called incipient spirals. In any 
event, should such "departures'* occur at rela- 
tively low altitudes, say traffic pattern 
height, it becomes difficult for the pilot to 
make adequate control adjustments. So-called 
"stall/spin" accidents are one of the greatest 
causes of fatalities in light aircraft. Engine 
failures on take-off, for example, are often 
followed by aircraft gyrations which inadver- 
tently precipitate severe wing lift loss. The 
propagation of the spanwise load changes on 
the stalling wing(s) is significant in the 
ensuing dynamics . 

This paper describes research which is 
to provide a logical basis for, and to demon- 
strate a method of developing relieving wing 
span-loadings which lead to a reduction in the 
post-stall departure forces. (Ultimately, of 
course, the research must include an analysis 


of the complete airplane dynamics problem to 
be valid.) 

BACKGROUND 

The parameters which directly influence 
the post stall handling qualities and depar- 
ture resistance of aircraft are listed below, 
from Ref. (1)*. 

1. Lift Curve Top, Cl versus a 

2. Rolling Moment, Ci 

3. Roll Damping, Cip 

4. Roll Control Power, Cl^^^ 

5. Aileron Adverse Yaw, Cn^ 5 ^ 

6. Directional Stability, Cjqg 

7. Yaw Damping, Cm^ 

8. Dihedral Effect, C^g 

9. Side Force, Cy, and Yawing 
Moment , C^ 

10. Yawing Moment Due to Rolling, C^p 

11. Pitching Moment, Cjq 

In general, the post-stall wing perfor- 
mance affects each of these parameters; some 

*Numbers in parentheses designate 
References at end of paper. 


ABSTRACT 


Results of recent wind tunnel tests and 
analyses are described which are aimed at helping 
to prevent spin entry after a stall by modifica- 
tions to the aircraft wing. Simple add-on devi- 
ces were utilized which modify the shape of the 
lift curve after stall, making it essentially 


flat to 30-40^ angle of attack. It is hypothe- 
sized that this will greatly reduce the danger 
of spin entry. Further research will be aimed 
at laying a firmer theoretical and analytical 
basis for the approach and at Incorporating the 
complete full scale airplane and its dynamics 
into the problem. 
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to a greater extent than the others. However, 
if the post-stall change in wing span loading 
could be designed to be gradual and nearly 
symmetrical about the aircraft roll and yaw 
axes, very gentle adverse accelerations would 
occur. Also, if adequate stall warning and 
good post-stall control could be sustained the 
research objectives would be met. It may thus 
be argued that the first step in solving the 
departure problem is to cause the post-stall 
spanwise wing unloading to be symmetrical and 
gradual, even under the adverse condition of 
a severe, tip-to-tip angle-of-attack variation 
(this has been the primary object of the research 
conducted thus far). For example, if one wing 
were stalled at a = 17° (average), and the other 
at a = 25°, no sudden roll accelerations should 
occur. This leads to the requirement for the 
so-called "flat top" lift curve. Now, if the 
drag difference between the two stalled wings 
is small and effective roll control without 
adverse yaw is available, many of the departure- 
inducing derivatives are significantly reduced. 

Though the shape of the three-dimensional 
wing lift curve proposed does not positively 
preclude fully developed spins (flat spins in 
particular) its development will certainly have 
a preventive influence on the initial post-stall 
acceleration rates and the pilot’s lateral, 
directional, and longitudinal control authority. 

Several practical constraints are always 
imposed on the solutions and thus the direction 
of research conducted for generalized light air- 
craft. These include limitations of weight, 
cost, simplicity, durability, reliability, main- 
tenance, manufacturing, performance and customer 
acceptance. Though flight safety is the ulti- 
mate goal of this research, only those solutions 
which satisfactorily meet the constraints are 
reported in this paper. 

Finally, the material presented has not 
been verified by full scale testing as of this 
writing. Its significance, however, v^as felt 
to be sufficient to warrant publication at this 
time . 

TECHNICAL DISCUSSION 

A combined theoretical and experimental 
aerodynamic research effort has been conducted. 

The theoretical analysis was based on the work 
reported in Refs. (2) through (4). The experi- 
mental work was performed in The University of 
Michigan and NASA Ames low-speed wind tunnels. 

It has been partially documented in Refs. (5) 
and (6). 

THEORETICAL APPROACH — Nonlinear lifting 
line theory was selected as the basis of the 
theoretical approach used in the research 
reported here. The primary constraint of the 
theory was that the post-stall nonlinear aero- 
dynamic characteristics of the elements (chord- 
wise truncated segments of the wing span) were 
known. More sophisticated nonlinear analyses 


have been reported, for example, Ref. (7). 
However, the extent of the development of Ref. 
(4) made the more simple approach most attrac- 
tive. 

In essence, the basis of the lifting line 
theory is as follows: If the wing may be 

represented by a set of bound vortex elements, 
having strengths which collectively provide the 
proper spanwise loading, then the local lift 
coefficient, c£(y,ag), chord, c(y), and span b , 
produce a spanwise downwash along the lifting 
line. This downwash is given by; 

V 

-b/2 

where is the local effective angle-of-attack, 
Voo is the free-stream velocity and y and y' are 
the spanwise coordinates. The local effective 
angle-of-attack becomes; 

a^(y) = a(y) - tan"^ 

00 

. ct(y) - ^ 

OO 

and the local lift coefficient becomes 

= f(%) 

Similarly the remaining nonlinear load coeffi- 
cients may be calculated once the effective 
local angle of attack is known. 

Figure 1 shows the types of section stall 
used in the analysis. Practically, the trailing- 
edge stall is accompanied by gradual turbulent 
separation from the rear of the airfoil, leading- 
edge stall is initiated by a short bubble from 
the leading edge and thin airfoil stall with a 
long bubble. The limitations of the theory 
were clearly demonstrated when the post-stall 
lift characteristics were varied. A doctoral 
dissertation issued by The University of 
Michigan (8) provides the proof of the existence 
of a mathematical singularity which arises as a 
consequence of the post-stall negative lift- 



SeCTION ANGLE OF ATTACK, 0 
Fig.l - Types of stall 



3 


curve slope. Qualitative post-stall results 
were obtained for the lifting line theory using 
a time dependent solution. In addition, the 
existence of lift hysteresis was indicated. 

Such effects are well known in helicopter dy- 
namics (9) . 

The results of the theoretical study may 
be summarized as follows: 

1. The method of Refs. (2)-(4) is inade- 
quate for general post-stall calculations. 

2. Qualitative, 3-D results have been 
obtained using a ’’transient solution” approach. 

3. Hysteresis in Cl vs a does exist 
mathematically . 

4. A nicely rounded (pure T.E. stall) sec- 
tion lift curve does not preclude rapid post- 
stall lift loss , 

5. A more clear understanding of wing 
stall propagation is becoming available. 

6. Stable asymmetrical lift solutions 
exist beyond stall. 

Additional theoretical studies were made 
using the linear vortex-lattice method to deter- 
mine the influence of the location of the 
stalled portion of the wing. Both span-wise 
and chordwise separation were simulated by 
allowing lift loss at the local ’’stalled” 
region. It was seen that local leading edge 
type separation produced a much stronger in- 
duced angle-of-attack distribution than trail- 
ing edge separation (6). This study led to the 
logic which ultimately guided the experimental 
program. Figure 2 shows how the trailing edge 
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Fig. 2 - The stalling "strip" model 


and leading edge separation were mixed to 
localize the wing stall propagation. In essence, 
the premature leading edge stall over a limited 
spanwise region unloaded the remainder of the 
wing by decreasing the adjacent geometric 
angle of attack. This was provided by the 
local induced angles-of-attack. Further, the 


stability studies of the non-linear lifting 
line theory implied the need for a local post- 
stall positive lift curve slope (following the 
rapid lift loss initiated at stall) to limit 
the induced total wing lift loss. Figure 3 




WHERE a = a, + Aa, 

'x 's \ ' 

DUE TO STALL 
(UNSTABLE) 

STALL a, 

Fig. 3 - The strip theory stalling thesis (strips 
1 and 2 from Fig. 2) 


demonstrates this argument. As 0.2 becomes 
unstable, it continues to diverge in effective 
angle-of-attack until it reaches a positive lift 
curve slope. Simultaneously, cxl is diminished 
or maintained until the entire wing reaches a 
stable point. Thus, if the stalled region can 
be localized (and stabilized), higher geometric 
angles-of-attack of the wing can be achieved 
with a reduced wing lift loss rate, resulting 
in approximately constant Cl to a high a. 

Lift hysteresis is also a consequence of 
the spanwise induced angle-of-attack instabil- 
ity. Figure 4 demonstrates how this can occur. 
After stall, the strip effective angle-of-attack 
continues to increase until a stable point is 
reached, well up on the post-stall positive lift 
slope. Then, as the geometric angle-of-attack 
is reduced, stability persists until the insta- 
bility re-appears , though at a lower angle-of- 
attack than at the previous stable point. Again, 
this is due to the induced flow interactions 
between all spanwise elements of the wing. 

When the reattachment occurs, the geometric 
angle-of-attack can be below that of the initial 
stall. The inset ’’wing Cl” curve in Fig. 4 
shows the effect on the entire wing. 

EXPERIMENTAL STUDIES - The experimental 
program was closely aligned with the theoreti- 
cal studies. The main objective was to find a 
means for containing the separating region so 
that the flat-top lift curve could be produced. 
The approach taken was as follows : 

1. Tuft studies for post-stall comparison 
of the wind tunnel model with a full scale flight 
test wing, 

2. Flow visualization of upper surface 
flow to allow a qualitative understanding of 
the separation mechanics, 

3. Investigation of devices to localize 
the separated flow, 

4. Study wing hysteresis, 
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WING a 

Fig. 4 - Stalling hysteresis illustration 


5. Investigate the effect of mixing lead- 
ing edge and trailing edge separation on a wing, 

6. Develop a flat top lift curve, and 

7. Improve performance. 

The requirements of the research program 
as of this writing limited the studies to 
simple rectangular wings. In addition, retro- 
fit devices were first studied because of their 
timeliness as defined by minimal manufacturing 
changes . 

Two slightly twisted (1.5°), aspect ratio 
seven wings with an NACA 652“^15 airfoil were 
obtained from NASA Langley. One was modified 
for the University of Michigan’s 5’ x 7’ wind 
tunnel installation and the other for the 
NASA Ames 7’ x 10’ wind tunnel. The wings 
were tufted to study the nature of the flow 
separation. At the tufts exhibited 

moderate spanwise flow at the trailing edge 
(Fig. 5). As the angle-of-attack was increased 
slightly, leading edge (large bubble) separa- 
tion occurred. The induced effects caused an 
upper surface stagnation vortex set to be shed. 
The strength of the cell, illustrated in Fig. 6, 
implies the level of lift loss. Also, the weak- 
ness of the inboard induced counter-rotating 
cell set implies the absence of upper surface 
lift in that area. It was interesting to note 
the dynamics of the upper surface flow during 
the incipient stall. An instability was present 
which changed the tuft patterns alternately 
(and implicitly, the lift) between the pattern 
shown and one with the vortex cell moved fur- 
ther outboard. This instability caused signifi- 
cant model motions as well as changes in tunnel 
flow conditions. It was noted in the lift read- 
out that more of a bi-stability as opposed to 
an instability was present. When the model 
angle-of-attack was increased the violent nature 
of this instability disappeared. Finally as 
the complete upper surface was stalled, though 
turbulence persisted, the wide ranging lift 
fluctuations disappeared. 

Next, a series of photographs was taken to 
demonstrate the flow preference with increasing 
versus decreasing angles-of-attack. Figure 7 
depicts the results. Note the inset sketch 


a = 17® 



FLOW AT CtMAx (STABLE) 


a = 21’* 
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INCIPIENT STALL (UNSTABLE) 


a = 25“ 



POST STALL (STABLE) 

Fig. 5 - Stall progression, basic wing 




STAGNATION VORTEX 



Fig. 6 - Example of a separation vortex 


showing model attitude and wind direction. It 
became clear that with increasing angles-of- 
attack, at a moderate rate of approximately 
2/10 degrees per second, the lift chose the 
higher value. Read upper photos from left to 
right. As the angles-of-attack were reduced 
from 27 to 14 degrees, from right to left in 
the photos, a hysteresis was demonstrated. 
Between 18 and 22 degrees (stall being at 
approximately 17°, uncorrected) unique flow 
patterns developed. 

A series of wind tunnel force measurement 
studies was next made to evaluate various types 
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Fig. 7 - Demonstration of flow hysteresis 


of separation stabilization devices. These 
included : 

1. Various forms of vortex generators, 

2, Fences, both streamwise and canted, 

3. Flexible (floating) trailing edge de- 
vices, and 

4, Various forms of leading edge vortex 
generators and stall strips. 

The results were that all of these devices, 
when applied to a wing prone to leading edge 
separation, produced results which were far 
too subtle to provide improvements at high 
angles-of-attack. 

Next, devices which would have a much 
stronger influence on the character of the flow 
separation were investigated. These included: 

1. Leading edge slats, and 

2. Drooped leading edges. 

It is felt that a more superior method than 
either of the above is now available through 
Ref. (10). This approach uses localized thick- 
ening of the upper surface toward the leading 
edge. However, as will be shown, the two devices 
listed above were sufficient to demonstrate the 
flat-top lift curve desired. Later studies will 
include the approach of Ref. (10) as it provides 
a generalized flexibility, can produce better 
performance and blends well with the manufactur- 
ing requirements. 

Figure 8 shows the nomenclature and des- 
cribes the stripwise quartering of the leading 
edge devices used in the experimentation. Vari- 
ous combinations of leading edge device loca- 
tions were attempted as a means of obtaining 
localized trailing edge separation. Tuft studies 
were used in the beginning and force measurements 
later. Also during the NASA Ames wind tunnel 
tests, one of the leading edge segments was split 
into two 12 1/2% span lengths (A & B) for further 
refinement . 

First, slats were used. The results showed 
significant promise as the stall could be tai- 
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"'L.E. SEPARATION STRIP 


Fig. 8 - Wind tunnel model layout 


lored to cause a loss of only 5% of the lift and 
quickly resume a positive dCL/dc^ slope. This 
could be sustained to angles-of-attack in excess 
of 40®. The best location for the unslatted 
portion was position 2 as shown in Fig. 8. More 
details of the study are given in Ref. (6). 

Next, drooped leading edge segments pattern- 
ed after the GAW-1 airfoil (Fig. 9) were used. 
Since the upper surfaces of both airfoils are so 
similar, this approach promised good cruise per- 
formance and trailing edge separation in stall. 
Wind tunnel tests were made at NASA Ames to 
determine the best spanwise location and width 
of the undrooped section. The results showed 
location 2B(Fig. 8) to be the best position to 
remove the 12 1/2% drooped segment. The tuft 
photos of Fig. 10 show the significance of 
bounding the leading edge separation with trail- 
ing edge separation. Figure 11 further shows 
the stability of the flow pattern up to 40® 
angle-of-attack. Comparison of the lift and 
drag data between the plain wing and configura- 
tion 2B is shown in Fig. 12. The hysteresis of 
the plain wing was obtained by averaging the 
data of several runs. Note the significant 
reduction in the pro-spin qualities of both lift 
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consequently, no lift carryover across the fuse- 
lage (with an attached region at midspan on each 
wing and generally, a separated region at the 
tip). These phenomena can be observed by care- 
ful study of the post-stall tuft photos for the 
652-415 wing with and without the "2B configura- 
tion" L.E. devices. 

The final wind tunnel tests of this series 
was for the evaluation of the stall character- 
istics of the 17%, GAW-1 wing, rectangular plan- 
form, aspect ratio seven. The model was obtained 
from NASA, Langley and was Identical in construc- 
tion to the previously tested wing. The Reynolds 


WING WITH 2-B MODIFICATION 



CLEAN 652-415 WING 



Fig. 10 - Comparison of wing flows at ot = 20® 


and drag using the modified wing. Also, the 
flap top lift curve has been demonstrated be- 
yond 30® angle of attack. Equally as signifi- 
cant is the maintenance of lift effectiveness 
with aileron deflection. Figure 13 shows no 
severe lift loss and thus implies aileron effec- 
tiveness beyond 40® angle-of-attack. The values 
of 15® down and 30® up aileron chosen are 
typical of fleet aircraft. 

A simplistic explanation of the reason for 
the success of the "2B configuration" in produc- 
ing a flattop lift curve is as follows: Moving 

the "L.E. stall strip" location outboard from 
the root allows a lift carryover to survive 
across the fuselage and produces three lifting 
cells across the span at high a’s as opposed 
to the two normally achieved with the "stall 
first at the root" philosophy, i.e. the wing, 
at post stall a’s, is transformed into three 
low aspect ratio surfaces: one across the 

fuselage and on the wing roots, and one on each 
tip section (with a separated area in between) . 
In contrast, the "conventional" stall pattern 
results in separation in the wing roots, and 


number was 2 x 10^ in the Ames 7* x 10* tunnel 
as before. Though the value of for the 

wing was somewhat higher than the ^2-415, the 
post-stall lift curve was continuously diminish- 
ing and without the upturn at 30® angle-of-attack. 

The hysteresis levels appeared to be much lower 

as was the average rate of lift loss with angle- 

of-attack. Several types of stall strips were 

used to induce leading edge separation in the 

same position, 2B. The results are shown in 

Fig. 14 for comparison with the 652-415 leading 

edge. An improvement in the flatness of the 

GAW-1 lift curve was attained with the stall 

strips though it was accompanied with a severe k 

loss in This type of behavior had been 

expected as a result of the analytical studies. 

CONCLUSIONS AND RECOMMENDATIONS 

The conclusions, based on the work to date 
may be summarized as follows: 

1. The non-linear lifting line theory as 
classically presented is not generally useful 
for post-stall studies. 







7 



a=30® a=4o° 


Fig. 11 - "Configuration 2B" showing separation 
localization at high angles-of-attack 



Fig, 13 - Effect of aileron on "configuration 2B" 
lift curve 

the untreated wing, 

5. The pro-spin hysteresis effect has been 
diminished by proper mixing of leading and 
trailing edge stall. 

6. The flat-top lift curve is accompanied 
by good lift and drag characteristics with post- 



a, deg 


Fig. 12 - Demonstration of "flat top" lift curve 


2. Lifting line solutions have been ex- 
tended beyond stall using a transient solution, 

3. Hysteresis exists in lift and drag 
with angle-of-attack under the conditions 
observed. 

4. An experimental method, based on the 
logic of an analytical model, has been developed 
for obtaining a *’flat-top** lift curve with no 
loss in CLmax or measurable Cj) increase over 


stall aileron deflection. 

7. A nicely rounded two-dimensional lift 
curve, accomplished with trailing edge separa- 
tion, does not preclude rapid three-dimensional 
post-stall lift loss if the full span section 
is constant (6) . 

8, Mathematically, asymmetrical full span 
stable lift distributions exist beyond stall 
for simulated steady unaccelerated flight (6). 
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a. deg 

Fig. 14 - Comparison with "GAW-1 wing" lift curves 

It is anticipated that the research will 
be continued with the following objectives; 

1. Develop a digital computer solution to 
predict the general post-stall wing load distri- 
bution for research and design purposes. 

2. Extend the theoretical and experimental 
investigations to include wing geometric charac- 
teristics such as taper, twist, flaps, and sweep. 

3. Apply the new concepts and theoretical 
modeling capability to the evaluation of 
incipient spin handling qualities. 

4. Demonstrate the spin resistant quali- 
ties of the flat top lift curve by full scale 
tests . 

The recommendations listed above are the 
goals of ongoing research being jointly 
approached by NASA Ames and The University 
of Michigan under NASA sponsorship. 
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